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An  experlasnt  Is  described  by  vhlch  the  ephemeris 
of  a  near-synchronous  satellite,  designated  by  I967-66F,  has 
been  determined  from  passive  range  observations.  More  pre¬ 
cisely,  the  data  consists  of  the  measured  times  of  reception 
at  ground  tracking  stations  of  electromagnetic  signals  which 
are  radiated  from  the  satellite  at  time  Intervals  regulated 
by  an  ultrastable  free. running  oscillator.  The  ground 
tracking  stations  mofcer  four  and  are  located  at  Howard 
County,  Maryland  (U.8.A.),  Ml  saws  (Japan),  Smlthfleld 
(Australia)  and  San  Jose  doe  Campos  (Brazil).  An  ephemeris 
and  three  satellite  clock  paramo  tore  are  estimated  by  the 
procedure  of  differential  correction  In  the  method  of  least 
squares.  A  comparison  of  the  ephemeris  to  one  obtained  from 
the  Doppler  tracking  method  Is  presented.  This  Indicates 
that  the  accuracy  of  the  ephemeris  is  better  than  4.9  mllll- 
radians  rm^  vhlch  was  the  goal  of  the  orbit  determination 
effort. 
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INTRODUCTION 


AWUCO  PMVtICt  UaOHATOMV 


This  docunent  reports  on  an  experiment  to  determine 
the  efheaerldee  of  a  near-synchronous  satellite  from  obser¬ 
vations  which  are  known  as  passive  range  measurements.  A 
passive  range  tracking  system  consists  of  measurements  taken 
at  tracking  sites  of  the  times  of  arrival  of  electromagnetic 
signals  transmitted  at  regular  intervals  from  the  spacecraft. 
The  time  Interval  between  the  Instant  the  signal  leaves  the 
spacecraft  and  when  it  arrives  at  a  tracking  site  is  a 
measure  of  the  range  to  the  spacecraft.  This  system  differs 
from  the  usual  radar  range  measurement  system  in  that  it  is 
open  loop)  the  ground  stations  are  radio  passive.  The  typical 
radar  range  measurement  system  is  closed  loop  in  which  the 
ground  station  transmits  a  signal  which  is  reflected  (either 
actively  or  passively)  from  the  spacecraft.  Inherent  to  the 
passive  range  technique  Is  the  necessity  to  maintain  time 
synchronization  between  the  clocks  at  each  of  the  tracking 
sites  and  the  spacecraft.  It  Is  possible  to  maintain  syn¬ 
chronization  of  the  station  clocks  by  several  means.  These 
are,  l)  to  hand  carry  portable  clocks  from  station  to  station, 
2)  to  utilize  W  radio  transmission  from  timing  reference 
stations  such  as  the  WV  station  at  the  National  Bureau  of 
Standards  and  3)  to  use  the  timing  transmissions  from  the 
Navy  Navigation  Satellite  System.  The  calibration  of  the 
satellite  clock  Is  distinct  from  the  station  clock  calibration 
because  the  satellite  Is  Inaccessible  and  presumably  at  an 
unknown  location.  Consequently,  It  Is  necessary  to  estimate 
the  satellite  clock  parameters  as  part  of  the  orbit  determin¬ 
ation  procedure. 

The  DODGE  (Department  of  Defense  Gravity  Experiment) 
satellite  program  vas’de signed  primarily  to"demonstrate 
gravity-gradient  stabilization  near  synchronous  altitudes  and 
to  determine  the  adequacy  of  theoretical  analyses  by  correla¬ 
tion  between  digital  simulation  results  and  experimental 
attitude  data  [1].*  A  complete  description  of  the  mission  as 
well  as  the  satellite  and  ground  Instrumentation  Is  given  In 
12]. 


Brackets  denote  References  at  end  of  paper. 
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The  nominal  prelaunch  orbit  characteristics  were 
an  altitude  of  16,200  nautical  miles  and  an  Inclination  of 
7  degrees.  This  would  imply  an  orbital  period  of  approxi¬ 
mately  22.166  hours  and  a  period  relative  to  an  earth- fixed 
meridian  of  approximately  13.2  days.  At  a  latitude  of  39 
degrees,  which  is  the  location  of  the  satellite  command 
facility  at  the  Applied  Jhyslcs  Laboratory,  the  spacecraft 
would  be  In  constant  view  for  about  5.6  days,  reappearing 
after  7*6  days. 

The  accuracy  requirement  for  the  ephemerides  was 
arrived  at  after  considering  several  aspects  of  the  satellite 
mission.  A  posteriori  ephemerides  are  necessary  for  the 
stabilization  studies  and  the  reduction  of  telemetry  and 
video  data  especially  that  from  which  the  attitude  of  the 
spacecraft  Is  to  be  determined.  A  maximum  a  posteriori 
position  error  of  100  nautical  miles  (2*. 6  nllllradlans)  was 
established  as  a  realistic  goal  for  the  DODGE  satellite.  By 
present  day  technology  this  Is  not  a  demanding  criterion, 
but  the  experiment  was  to  examine  the  feasibility  of  the 
method  with  the  Intention  of  Improvements  on  subsequent 
spacecraft. 


-  2  - 


ARRtlfD  PHYSICS  LABORATORY 

lw«  Mat  MHHMH 


2.  TRACKING  METHOOOLOGY 


Wie  techniques  most  frequently  used  to  obtain 
satellite  tracking  data  can  be  grouped  Into  three  categories: 
angle,  range-rate  and  range.  There  are  several  different 
methods  by  which  data  of  each  type  can  be  obtained.  For  low 
to  intermediate  altitude  satellites,  data  from  any  one  of 
these  techniques  would  be  sufficient  to  establish  an  orbit. 
For  example,  the  Applied  Ihysics  Laboratory  has  developed 
the  Navy  Navigation  Satellite  System.  Here  the  measurement 
of  the  Doppler  effect  on  a  VHF  carrier  (range-rate)  pro¬ 
vides  the  data  from  which  the  satellite  epheraeris  is  deter¬ 
mined.  The  precision  of  the  a  posteriori  ephenerides  of 
these  navigation  satellites,  which  are  at  a  nominal  altitude 
of  600  nautical  miles,  is  better  than  20  meters  [3]. 

While  the  Doppler  method  works  well  for  low  alti¬ 
tude  satellites,  degradation  is  to  be  expected  for  high  alti¬ 
tude  satellites  where  the  magnitude  of  the  Doppler  effect  is 
reduced  and  the  jasurement  errors  remain  the  same.  In 
addition,  for  a  truly  geostationary  satellite,  there  would 
be  no  measurable  Doppler  effect  and  the  longitude  (along- 
track  position)  of  the  satellite  would  be  indeterminate. 

Because  of  the  existence  of  the  Navy  Doppler  track¬ 
ing  stations  and  the  expertise  realized  in  the  development 
of  the  Doppler  tracking  technique,  studies  were  first  under¬ 
taken  to  establish  the  merit  of  Doppler  tracking  of  the  DODGE 
spacecraft  (L),  These  studies  shewed  that  while  it  might  be 
possible  to  achieve  the  required  accuracy  of  100  nautical 
miles,  the  confidence  in  doing  so  was  not  high.  Conse¬ 
quently,  an  auxiliary  mean3  for  tracking,  designated  the 
DGDGE  Time  Recovery  System  (DTRS),  was  developed  (5). 

Basically,  the  DTRS  provides  satellite  instrumenta¬ 
tion  for  periodic  transmission  of  a  "time  pulse",  a  discrete 
identifiable  modulation  of  the  dcwn-llnk  carrier.  The  control 
of  this  transmission  is  provided  by  the  satellite’s  automatic 
programmer  which,  in  turn,  uses  the  satellite’s  ultrastatle 
oscillator  as  a  basic  time  and  frequency  reference.  Appro¬ 
priate  equipeaent  was  designed  and  installed  at  four  existing 
TRANET  stations  to  receive  these  special  transmissions. 
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3.1  Satellite  and  Station  Instrumentation 


Fig.  3.1  is  a  block  diagram  of  the  satellite 
Instrumentation  used  in  generation  and  transmission  of  the 

time  pulses. 

Hie  time  pulses  appear  as  a  60  degree  phase  modu¬ 
lation  of  the  240  MHz  downlink  carrier.  Specifically,  the 
time  pulse  waveform  for  any  individual  transmission  consists 
of  128  cycles  of  a  square  wave  of  frequency  325.5  Hz.  Hie 
first  126  cycles  have  one  phase.  After  the  126th  cycle, 
there  is  a  180  degree  phase  shift  followed  by  two  more  cycles. 
This  phase  shift  is  the  signal  feature  which  is  detected  and 
measured  by  the  ground  station  equipment.  This  waveform  is 
shewn  in  the  top  half  of  Fig.  3.2. 

The  timing  mark  generator  shown  in  Fig.  3.1  trans¬ 
mits  a  128  cycle  burst  of  393  milliseconds  duration.  This 
sequence  is  repeated  every  3.2  seconds.  The  automatic  pro¬ 
grammer  sequences  these  bursts  into  the  total  schedule  of 
dcwnlink  transmissions  by  controlling  the  modulation  index 
control  unit.  The  programmer  starts  automatically  once  per 
hour  and  runs  the  satellite  for  ten  minutes.  It  will  then 
shut  off  unless  commanded  otherwise.  During  this  10  minute 
period,  three  time  bursts  will  be  sequenced  into  the  trans¬ 
mission.  These  occur  at  the  end  of  consecutive  200  second 
intervals,  as  shown  in  the  bottom  half  of  Fig.  3.2.  There 
are  thus,  under  normal  automatic  operation,  three  discrete 
pulses  per  ho'ir  or  72  data  points  per  day. 

Hie  automatic  programmer  and  the  timing  mark  gen¬ 
erator  obtain  their  frequencies  by  dividing  down  the  signal 
from  the  stable  5  MHz  oscillator.  Hius,  the  time  interval 
between  consecutive  time  pulses  is  tied  directly  to  this 
oscillator  as  a  reference.  For  this  reason,  the  oscillator 
must  have  been  turned  on  and  allowed  to  stabilize  before 
useful  data  can  be  obtained. 

Each  tracking  station  participating  in  the  time 
pulse  tracking  network  is  equipped  with  a  unit  capable  of 
acquiring  the  time  burst  waveform  and  detecting  the  phase 
reversal.  This  detection  triggers  a  reading  of  the  station 
ultrastable  clock.  This  reading  is  the  basic  measurement. 

frictfiig  pi|>  Milk 
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Fig.  3.1  DOOGE  SATELLITE  TIME  PULSE  SYSTEM 


6 


t m  mm  no— 1«§  iww 
Am.)KO  PHYSIC*  LASOMATOHY 
•urn  iMMa  Munaw 


Corrected  for  measured  delay  in  the  receiving  equipment  and 
for  epoch  and  drift  in  the  station  clock,  the  reading  becomes 
a  time  pulse  data  point.  The  DTRS  ground  system  tests  are 
described  in  detail  in  [6]. 

Four  stations  of  the  TRANET  system  have  been 
equipped  with  these  DTRS  units:  station  008j  j»an  Jose  dos 
Campos,  Brazil;  station  013,  Misava,  Japan;  station  112, 
Smithfield,  Australia;  and  station  §11,  Howard  County, 
Maryland*  The  units  are  able  to  record  the  time  of  recogni- 
trion  of  the  time  burst  from  the  satellite  with  a  one  micro¬ 
second  resolution. 

An  obvious  basic  requirement  is  that  the  station 
clocks  in  this  network  be  synchronized  to  high  precision. 

Uiis  means  that  an  epoch  must  be  established  at  each  station 
with  respect  to  a  common  reference.  Similarly  a  common 
clock  rate  reference  must  be  used. 

The  epoch  is  best  established  by  periodic  visits 
to  the  station  with  a  portable  atcsnic  clock  which  yields 
time  transfers  accurate  to  a  few  microseconds.  A  description 
of  this  process  and  results  of  visits  to  the  stations  are 
discussed  on  pp.  43-50  of  [7].  The  frequency  or  clock  rate 
reference  is  obtained  by  station  tracking  of  standard  VLF 
transmissions.  For  the  time  pulse  experiment,  the  Pacific 
stations  at  Smithfield  and  Mi6awa  used  the  13.6  KHz  trans¬ 
missions  from  Haiku,  Hawaii;  this  is  an  OMEGA  station.  The 
station  at  San  Jose  dos  Campos,  used  the  12  KHz  OMEGA  signal 
from  Trinidad,  while  the  Howard  County  clock  is  maintained 
in  synchronism  with  the  U.S.  Naval  Observatory,  A  discussion 
of  this  frequency  tracking  and  associated  instrumentation  at 
a  station  is  presented  in  Section  2,  Chap.  9  of  [5]. 

3.2  System  Error  Sources 

There  are  a  number  of  errors  in  the  time  pulse 
tracking  process  described  above  which  will  cause  errors  in 
the  orbit  determined  from  such  data.  In  this  subsection, 
these  are  summarized,  leaving  the  mathematical  representation 
for  the  next  subsection. 

First,  there  is  the  stability  of  the  satellite 
instrumentation,  particularly  its  5  Mz  ultrastable  oscilla¬ 
tor.  It  is  also  possible  to  interrupt  the  regular  sequence 
of  time  pulses  by  commands  which  reset  the  automatic  pro¬ 
grammer  or  switch  between  the  pair  of  oscillators.  These 
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are  discrete  events  that  are  easily  recognizable  In  the  data 
through  the  large  discrepancies  in  programmed  Interpulse 
Intervals  which  they  cause.  However,  they  must  be  guarded 
against  because  they  break  up  the  span  of  data  that  can  be 
employed  In  tracking. 

The  ultrastable  oscillator  units  used  on  DODGE  are 
of  the  sane  quality  as  those  used  in  the  Navy  Navigation 
Satellite  System.  They  can  thus  be  expected  to  maintain 
short-term  drift  rates  of  a  few  parts  in  10^  and  long-term 
drift  rates  of  1  or  2  parts  in  10^  per  day.  Normally,  one 
would  not  expect  such  drift  rates  to  be  constant  over  long 
periods  of  time,  i.e.,  days  or  weeks.  Rather,  one  would 
expect  them  to  have  something  of  a  "random-walk"  behavior, 
changing  in  direction  and  magnitude  at  random  and  unpre¬ 
dictable  times. 

A  drift  rate  of  2  parts  in  10*°,  maintained  over 
a  one- day  interval,  would  result  in  a  9  usee  error  in  pro¬ 
grammed  time  of  transmission  for  a  pulse.  Because  of  the 
random-walk  character  of  the  drift,  it  was  originally  hoped 
that  an  error  of  this  size  would  be  rather  infrequent.  How¬ 
ever,  as  will  be  discussed  later,  the  data  utilized  indicated 
that  the  satellite  oscillator  drift  was  apparently  nearly 
constant  over  the  whole  tracking  span.  Consequently,  it 
became  necessary  to  account  for  this  drift. 

The  ionosphere  introduces  variations  in  the  group 
path  velocity  of  the  VHF  carrier.  These  introduce  errors 
in  the  measured  times  of  arrival  of  the  pulses.  The  amount 
of  this  error  is  highly  dependent  on  the  profile  of  electron 
density  in  the  ionosphere,  and  also  on  the  elevation  of  the 
satellite  with  respect  to  the  station.  The  electron  density 
profile  has  a  high  degree  of  variability  on  a  diurnal,  a 
seasonal  and  a  geographical  basis. 

To  study  the  ionospheric  refraction  error,  several 
profiles  were  postulated,  representing  both  "ordinary"  and 
"extreme"  conditions.  The  equations  representing  the  transit- 
time  error  for  240  Miz  signals,  in  terms  of  each  such  profile, 
were  then  numerically  Integrated.  The  results  furnish  a 
spectrum  of  possible  ionospheric  errors.  These  studies  are 
reported  in  detail  in  [8].  Suffice  it  to  say  that  the  normal 
ionospheric  error  is  in  the  range  of  a  few  microseconds  and 
hence  was  neglected  in  the  orbit  determination  process. 
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The  remaining  system  errors  are  connected  with 
station  instrumentation  and  operation.  They  fall  into  three 
principal  areas: 

(i)  The  accuracy  with  which  remotely  separated 
clocks  can  be  synchronized  to  the  common 
UTC  standard. 

(ii)  The  accuracy  with  which  equipment  delays  can 
be  measured  and  controlled, 

(iii)  The  signal-to-noise  ratios  in  the  DTRS  units. 

These  classes  of  errors  are  discussed  fully  in  Section  2, 

Chap.  9  of  [5]*  That  discussion  is  summarized  here. 

The  r.f.  noise  error  (iii)  can  be  minimized  by 
employing  suitably  narrow  bandwidths  in  the  tracking  filter 
of  the  DTRS  unit.  The  phase  reversal  of  the  input  signal  is 
detected  by  comparing  the  input  and  output  of  the  filter. 

Due  to  the  narrow  filter  bandwidth,  the  output  signal  con¬ 
tinues  at  its  original  phase  immediately  after  the  phase 
reversal  and  hence  the  input  and  output  signals  differ  by 
180  degrees.  A  phase  reversal  detector  compares  the  filter 
input  and  output  signals  and  triggers  a  time  measurement  on 
the  next  filter  output  cycle.  In  this  way,  a  white  noise 
bandwidth  corresponding  to  the  filter  bandwidth  is  achieved. 

A  calculation  given  on  pp.  9-35  of  [5]  indicates  that,  using 
a  filter  installed  in  a  temperature  controlled  oven  with  a 
bandwidth  of  2.4  Hz,  the  r.f.  noise  error  should  have  an 
rms  of  at  most  a  few  microseconds  for  the  lowest  input  signal- 
to-noise  ratio. 

The  design  objective  for  the  synchronization  error 
(l)  was  set  at  20  usee  peak-to-peak,  i.e.,  a  bias  error  of 
no  more  than  10  usee  relative  to  UTC  at  any  station.  This 
objective  requires  careful  setting  of  an  epoch  by  means  of 
a  portable  clock,  as  well  as  consistent,  accurate  monitoring 
of  the  VU  frequency  reference.  A  time  and  frequency  control 
experiment  carried  out  between  a  TRANET  site,  station  103, 

Las  Cruces,  and  station  911  j  Howard  County,  indicated  that 
such  accuracies  could  be  obtained  in  the  continental  U.S. 
However,  it  is  less  certain  how  well  they  could  be  achieved 
on  a  world-wide  basis.  This  Is  discussed  on  page  50  of  [7]. 

There  are  additional  problems  when,  for  various 
reasons,  a  station  loses  an  epoch  which  has  been  established 
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by  the  portable  clock.  A  careful  examination  of  station 
records  for  frequency  tracking  is  then  needed  to  infer  the 
probable  epoch  error.  It  is  conjectured  that  epoch  cannot 
consistently  be  established  much  better  than  50  iisec  under 
these  conditions. 

The  stations  in  the  DQDOE  time  tracking  network 
reset  their  clocks  daily  during  tracking  periods.  Examina¬ 
tion  of  the  frequency  errors,  as  determined  by  comparison 
against  the  VU1  references,  consistently  indicates  that 
accuracies  better  than  1  part  in  10^-°  have  been  maintained. 
Thus,  it  seems  likely  that,  from  the  viewpoint  of  frequency 
tracking  alone,  clock  errors  should  be  below  the  10  usee 
objective.  The  total  epoch  error  may,  however,  be  several 
times  this  amount. 

Hie  remaining  station  error  (ii)  results  from 
calibration  of  signal  delay  in  the  DTRS  unit,  nils  is 
discussed  in  [5]  and  also,  in  more  detail,  in  [9].  Tests 
were  made  to  determine  equipment  delay  times  as  a  function 
of  r.f.  signal  level,  receiver  loop  bandwidth,  receiver 
tuning,  AC  line  voltage  variations,  and  temperature  varia¬ 
tions.  The  results  in  [9]  indicate  that  the  dispersion  in 
successive  delay  measurements  is  about  6  usee,  nils  repre¬ 
sents  variability  in  delay  with  respect  to  a  dally  mean 
value.  In  addition,  there  is  some  long-term  drift  in  delay 
values.  Except  for  the  station  at  Mlsawa,  delay  measurements 
indicate  long-term  drifts  of  less  than  one  microsecond  per 
day. 


Summarizing  this  section,  the  major  errors  in  the 
measurement  of  time  points  arise  from  satellite  oscillator 
bias  and,  occasionally,  from  consistent  oscillator  drift, 

To  a  lesser  extent,  error  is  Introduced  by  inaccuracies  in 
synchronizing  all  stations  to  a  common  time  base  and  by  lack 
of  precise  signal  delay  calibrations  at  the  various  stations. 

In  principle,  one  would  hope  to  compensate  for 
satellite  frequency  offset  and  drift  by  solving  for  these  as 
part  of  the  orbit  determination  process.  However,  it  must 
be  remembered  that,  for  a  given  tracking  station,  a  DODGE 
"pass"  Is  4-5  days  In  duration.  It  seems  unlikely  that 
typically  the  oscillator  would  have  a  constant  drift  rate 
over  so  extended  a  period.  Hence,  to  hypothesize  one  and 
attempt  to  solve  for  it  Involves  a  theoretical  error.  It 
Is  also  true  that  the  least  squares  normal  equations  used  In 
orbit  determination  become  quite  ill-conditioned  when  drift 
rate  la  included  as  a  parameter.  Nor  Is  the  amount  of  data 
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available  In  a  single  pass  from  a  single  station  ordinarily 
enough  to  give  even  a  good  determination  of  frequency  offset. 
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4.  ORBIT  DETERMINATION  PROGRAM 


4,1  Mathematical  Formulation 

We  consider  a  time  pulse  vhose  distinctive  feature , 
the  phase  reversal  of  Fig,  3.2,  is  emitted  from  the  satellite 
at  UTC  Instant  t .  .  Let  this  pulse  be  detected  by  station 
k  at  UTC  Instant1  .  This  Implies  that  all  necessary 
corrections  for  station  clock  error,  signal  delay,  etc., 
have  been  made.  Then,  since  refraction  effects  are  neglected. 


where 


pik  ■  i?(V  -  Wvi 


(k.l) 


and  r(t,)  Is  the  satellite  position  at  t,  and  ^.(t  . ) 
is  the  station  position  at  *  ' 

It  Is  clear  from  the  discussion  In  subsection  3.1 
that  there  are  two  kinds  of  intervals  between  any  two  con¬ 
secutive  time  pulses.  One  is  the  "small''  Interval  of  approx¬ 
imately  200  seconds  between  consecutive  bursts  within  a  given 
hourly  transmission.  The  other  Is  the  longer  Interval  (about 
3200  seconds)  between  the  last  pulse  of  one  hour  and  the 
first  pulse  of  the  next  hour.  Let  us  denote  the  exact  dura¬ 
tions  of  these  two  kinds  of  Intervals  by  Ajt  and  Agt  , 
respectively. 

Since  ve  do  not  know  the  exact  time-varying  history 
of  the  satellite  oscillator's  frequency  (assuming  again  that 
It  behaves  something  like  a  random  walk)  we  will  not  know 
A]t  or  Apt  .  However,  ve  can  form  estimates  of  them  and, 
perhaps,  at  least  partially  correct  these  estimates  by  means 
of  the  data.  Specifically,  let  our  best  estimate  of  the 
oscillator  frequency  In  Hz  be  denoted  f_  .  Since  there  is 
a  precise  known  number  of  oscillator  cycles  In  each  of  the 
two  kinds  of  Intervals,  the  current  best  estimates  of  the 
duration  of  these  Intervals  are 


13 


APPLIED  PHYSICS  LAaCMATOftV 


(Axt)e  =  (l, 006,632, 960)/fn 
(^t)e  .  (I5,986,196,473)/fn 


(4.2) 


The  numerators  are  the  cycle  counts  built  into  the  timing 
mark  generator  and  automatic  programmer.  The  subscript  "e" 
denotes  estimate.  Insofar  as  the  true  behavior  of  the  sat¬ 
ellite  frequency  over  the  tracking  span  can  be  characterized 
as  a  simple  frequency  offset  of  Af  Hz,  an  Improved  estimate 
would  be 


-  (41t)e  (1  -$£) 

"  \  <<*.3) 

45 1  -  (V>e(l-r> 

n  / 

Returning  to  the  basic  data  eq.  (4.1),  we  see  that 
we  can  replace  the  left-hand  side  with 

t4  -  t0  ♦  6t  +  k(i,l)  A*t  +  k(i,2)  A^t  (4.4) 


In  writing  this  equation,  ve  are  considering  that  data  for 
a  given  tracking  span  has  been  chronologically  arranged  and 
the  Initial  pulse  Is  estimated  to  have  been  transmitted  at 
UTC  Instant  tQ  .  There  will  be  an  epoch  error  In  this 
estimate,  of  course,  which  has  been  denoted  5t  in  eq.  (4.4). 
After  the  Initial  pulse,  all  subsequent  pulses  must  then  have 
been  transmitted  at  some  Integer  combinations  of  Ajt  and 
Agt  Intervals.  For  the  1-th  pulse  In  the  tracking  span  these 
Integers  have  been  denoted  k(l,l)  and  k(l,2)  .  They  are, 
of  course,  easily  determined  from  the  data  t..  since 

6t  «  Ajt  . 

Eqs.  (4.l)-(4,4)  are  the  theoretical  model  of  the 
time  pulse  data.  They  Involve  the  six  orbital  elements 
through  the  p~  and  the  two  "clock  parameters"  6t  and 
Af  ,  a  total  or  eight  parameters  to  be  determined  from  the 
data.  Simulation  Indicated  that  this  Is  the  maximum  one 
may  reliably  hope  to  determine  by  straightforward  least 
squares  procedures  using  tracking  spans  of  3*6  weeks  with 
each  station  averaging  60<  recovery  of  available  time  pulses 
(typical). 
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In  the  orbit  adjustment  program  discussed  in  sub¬ 
section  4.3  we  require  the  partial  derivatives  of  the  data 
t„  with  respect  to  the  8  parameters.  The  partials  with 
respect  to  Af  and  6t  are  correctly  determined  only  if 
we  remember  that  a  change  Af  to  f  does  two  things: 


(a)  It  changes  the  absolute  UTC  instants  of  all 

pulses  and  hence  also  the  satellite-station 
geometry  for  any  given  pulse,  i.e.,  . 

(b)  It  changes  A*t  and  A*t  . 


From  these  facts  we  find 


ik 


■STst) 


=  1  + 


ik 


(4.5) 


where  is  the  relative  range-rate  at  ti  .  Similarly 


*ik 


P41c  k(i,l)  Art  +  k(i,2)  A*t 

(1+-^)  - S - -  (4.6) 
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Let  K  denote  any  one  of  the  six  orbital  elements. 
Then  we  find  the  partial  by  the  usual  chain-rule  logic: 


*1* 

"ST  * 


l  *1* 

c  wr 
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(4.7) 


Here  denotes  the  (l  x  3)  row  vector  whose  elements 

are  the  partials  of  range  with  respect  to  Cartesian  satellite 
position,  all  evaluated  at  t^  ;  the  5X^/5K  is  the  (3  x  l) 
column  vector  whose  elements  are  the  partials  of  the  Cartesian 
satellite  position  components  (at  t.)  with  respect  to  the 
orbital  element  K  . 


4,2  Editing  the  Data 

Before  we  attempted  a  differential  correction  or 
least-squares  solution  for  the  8  parameters  discussed  in  the 
previous  subsection,  it  was  necessary  to  delete  obviously 
bad  data  points.  No  entirely  satisfactory  method  for  doing 
this  has  been  found.  The  method  which  finally  evolved  will 
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be  described  here  without  any  claim  that  it  is  in  any  sense 
"best" , 

Since  there  are  two  tracking  stations  in  the 
western  hemisphere  and  two  in  the  Pacific,  DODGE  was  nor¬ 
mally  visible  at  two  of  the  four  stations  at  any  one  time. 
Consequently,  it  was  possible  to  arrange  the  data  for  a 
given  tracking  span  into  "passes"  and,  for  each  pass,  to 
list  out  those  pulses  received  at  both  of  the  stations  then 
observing  the  satellite. 

This  procedure  involves  ignoring  a  significant 
amount  of  the  total  data  in  practice  because  there  are  many 
pulses  received  at  one  of  the  two  stations  but,  for  various 
reasons,  not  at  the  other.  However,  the  editing  operates 
only  on  those  pulses  received  at  both  stations. 

Using  our  best  current  estimate  of  the  orbit,  we 
compute  theoretical  arrival  times  for  the  pulses  at  each 
station  and  form  the  residuals.  We  then  subtract  the  resi¬ 
dual  for  one  station  from  the  residual  for  the  other.  This 
quantity  is  referred  to  as  "delta-delta".  Now,  although 
the  residuals  themselves  are  fairly  large  and  vary  with  time, 
the  delta-delta  quantities  effectively  subtract  out  any  time 
and  frequency  errors  involved  in  the  theoretical  data,  i.e., 
the  quantity  tj.  of  eq.  (4.4),  leaving  only  the  small 
timing  errors  due  to  satellite  orbit  errors.  However,  an 
error  in  the  experimental  data  that  occurs  at  one  of  the 
two  stations  will  appear  in  full  in  delta-delta.  This 
statement  of  course  ignores  the  possibility  of  compensating 
errors  occurring  at  the  two  stations, 

By  visual  scanning  of  the  delta-deltas  as  a  function 
of  time,  one  can  note  any  outliers  or  wild  points  and  manually 
delete  them  from  the  data.  Usually,  by  comparing  residuals 
from  each  station  before  and  after  the  point  in  question,  one 
can  assess  fairly  well  the  station  at  fault.  Although  not 
a  hard  and  fast  rule,  a  jump  of  100  microseconds  or  more  in 
delta-delta  has  been  used  as  a  basis  for  deletion. 

4,3  Orbit  Adjustment  Program 

The  final  result  of  the  editing  process  described 
in  subsection  4,2  above  is  a  tape  containing  that  part  of  the 
data  taken  during  the  tracking  span  which  is  "good  enough" 
to  furnish  an  orbit  correction.  To  process  this  tape  and 
derive  the  orbit  correction  is  the  function  of  the  orbit 
adjustment  program. 
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Like  most  orbit  correction  programs,  this  one  makes 
iterative  passes  over  the  data,  using  as  current  estimate  of 
the  8  parameters  the  values  derived  on  the  previous  iteration. 
The  initial  estimate  of  the  epoch  error  6t  is  zero.  For 
the  most  part,  the  initial  estimate  of  the  frequency  offset 
Af  has  been  obtained  by  making  one  pass  over  the  data  using 
Af  =  0  and  then  estimating  a  crude  value  of  Af  that  will 
approximately  null  out  the  predominant  linear  growth  of  the 
residuals  in  time. 

On  each  iteration  of  the  orbit  adjustment  program, 
processing  begins  by  converting  the  current  best  estimates 
of  Kepler  elements  at  the  epoch  to  an  internal  set  of  ele¬ 
ments  described  in  [U],  Numerical  integration  involving  the 
geopotential  forces  of  the  earth  complete  through  the  fourth 
order  and  degree,  the  perturbing  forces  of  the  sun  and  moon 
and  the  force  of  the  solar  radiation  is  used  to  generate  the 
ephemeris  from  which  the  "theoretical  values"  of  the  obser¬ 
vations  can  be  computed.  The  normal  equation  which  arises 
in  the  least-squares  method  is  solved  by  a  routine  discussed 
in  [10].  Because  the  problem  is  nonlinear,  it  is  necessary 
to  iterate  the  solution.  Breakout  occurred  after  a  fixed 
preset  number  of  iterations  so  that  the  parameters  could  be 
hand  adjusted,  if  necessary,  before  continuing  the  iteration 
process. 
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5.  DESCRIPTION  OF  THE  EXPERIMENT 


The  DODGE  spacecraft  was  successfully  launched 
into  orbit  on  1  July  1967  and  designated  as  I967-66F.  Hie 
orbital  characteristics  that  were  attained  were  close  to  the 
nominal  with  an  apogee  of  I8,l6l  nautical  miles,  a  perigee 
of  17,959  nautical  miles  and  an  inclination  of  7.0^7  degrees. 
Both  Doppler  and  time  pulse  tracking  measurements  were 
obtained  by  the  tracking  stations  indicating  that  the  hard¬ 
ware  in  the  spacecraft  and  at  the  stations  were  functioning. 
Because  the  Doppler  technique  had  been  designated  as  the 
primary  tracking  method,  attention  was  first  directed  toward 
it.  The  results  of  the  orbit  determination  procedure  utili¬ 
sing  the  Doppler  data  are  given  in  [;*]  where  it  is  demon¬ 
strated  that  the  accuracy  achieved  was  better  than  the  goal 
of  100  nautical  miles  rms.  Both  antenna  pointing  data  and 
measurements  of  the  time  of  passage  of  the  spacecraft  into 
or  out  of  the  penumbra  of  the  earth  were  utilized  to  obtain 
the  measure  of  accuracy. 

Study  of  the  time  pulse  observations  (i.e.,  the 
difference  between  the  time  of  reception  of  the  time  pulse 
and  the  time  of  transmission  from  a  nominal  satellite  clock, 
as  a  function  of  universal  time)  disclosed  a  phenomenon  in 
the  character  of  the  observations  that  had  not  been  expected. 
Since  the  satellite  moves  slowly  relative  to  a  point  fixed 
on  the  earth,  the  data  should  also  be  both  slowly  varying  and 
a  continuous  function  of  time,  During  tracking  periods  in 
1968  and  1969,  data  obtained  in  the  eastern  hemisphere  by 
the  stations  at  Misawa  and  Smithfield  did  indeed  have  this 
character.  However,  the  data  obtained  in  the  western  hemis¬ 
phere  at  both  Howard  County  and  San  Jose  dos  Campos  were  of 
a  distinctly  different  nature  exhibiting  Jump  discontinuities 
and  rapid  changes  in  amplitude.  Since  similar  anomalies 
appeared  in  the  observations  at  both  stations  at  the  same 
time  the  source  of  the  difficulty  had  to  be  at  the  satellite 
and  not  at  the  stations.  Further  study  uncovered  the  fact 
that  each  anomaly  was  correlated  in  time  with  commands  trans¬ 
mitted  to  the  spacecraft  from  the  command  facility  at  Howard 
County.  During  periods  in  which  there  were  no  spacecraft 
commands  transmitted,  the  observations  were  well-behaved. 

The  mechanism  by  which  the  performance  of  the  timing  mark 
generator  was  affected  by  commanding  could  not  be  uncovered. 
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Since  the  orbit  determination  procedure  utilizing 
the  Doppler  data  yielded  ephemerides  of  sufficient  accuracy, 
there  was  little  need  for  time  pulse  tracking  to  carry  out 
the  mission  of  the  spacecraft.  In  order  to  perform  a  defini¬ 
tive  experiment  utilizing  the  time  pulse  data,  it  was  nec¬ 
essary  to  wait  for  a  period  of  time  in  which  it  would  be 
possible  to  refrain  from  transmitting  commands  to  the  space¬ 
craft.  Such  an  opportunity  occurred  during  a  51-day  interval 
in  1970 >  from  day  17’*  to  day  225.  During  this  period  safe¬ 
guards  were  employed  to  ensure  that  the  spacecraft,  which  was 
in  the  automatic  programmer  mode,  would  not  be  accidentally 
commanded.  Because  the  maintenance  of  time  synchronization 
between  the  four  tracking  stations  is  imperative  for  meaning¬ 
ful  time  pulse  data,  a  concerted  effort  to  do  so  was  under¬ 
taken.  Estimates  of  the  rms  accuracy  of  the  a  posteriori 
calibrations  were:  5  ^sec  at  Howard  County,  20  ^see  at 
Smithfield,  50  nsec  at  Misawa  and  100  usee  at  S'  n  Jose  dos 
Campos.  The  calibration  of  the  clock  at  San  Jose  dos  Campos 
unfortunately  had  to  be  accomplished  without  the  benefit  of 
clock  transfers  that  were  utilized  at  the  other  sites. 
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6.  EFHEMERIS  DETERMINATION 


A  total  of  20  passes  of  data  Van  obtained  by  the 
four  tracking  station.;  during  the  51  day  experiment.  An  a 
posteriori  adjustment  was  made  to  the  data  utilizing  the 
best  estimates  of  the  offsets  of  the  station  clocks.  The 
data  obtained  at  all  the  stations  except  at  Howard  County 
when  plotted  as  a  function  of  time  were  continuous  and  slowly 
varying.  The  data  obtained  at  Howard  County  again  was  of 
poor  quality.  The  large  scatter  asso"4atel  with  this  data 
was  attributed  to  poor  signal  strength.  Since  the  site  at 
San  Jose  dos  Campos,  which  received  the  same  time  pulse 
signals  as  Howard  County,  obtained  the  signals  at  normal 
signal  strengths,  it  appeared  that  the  liff lenities  had  to 
be  in  the  receiving  equipment  at  Howard  County,  Several 
attempts  were  made  to  localize  the  difficulty  but  these 
proved  unsuccessful.  Subsequently,  the  observations  obtained 
at  Howard  County  were  not  used  so  that  the  total  population 
of  tracking  data  consisted  of  the  li*  passes  obtained  ut  Misuwa, 
Smithfield  and  San  Jose  dos  Campos, 

By  not  utilizing  the  Howard  County  data,  it  was 
necessary  to  suspend  the  editing  test  for  the  data  taken  at 
San  Jose  dos  Campos.  Instead,  only  obviously  bad  data  were 
deleted  by  "hand". 

Initial  attempts  to  fit  both  an  ephemerio  and  sat¬ 
ellite  oscillator  frequency  resulted  in  rather  poor  residuals. 
Study  of  the  behavior  of  the  residuals  suggested  that  there 
was  significant  drift  of  the  satellite  oscillator  frequency. 
Since  a  frequency  drift  parameter  was  not  one  of  the  original 
eight  fitting  parameters,  it  was  necessary  to  introduce^it 
as  the  ninth.  This  would  effectively  make  Ajt  and  A_t 
linear  functions  of  time.  To  minimize  the  magnitude  of*- 
the  modification  of  the  software,  the  adjustment  of  the  fre¬ 
quency  parameter  was  not  made  part  of  the  least- squares 
fitting  procedure.  Instead,  it  was  added  as  a  parameter 
which  could  be  adjusted  prior  to  each  iteration  of  the  least- 
square  process.  Sequentially  fitting  the  frequency  drift 
parameter  with  the  other  eight  parameters  resulted  in  mini¬ 
mizing  the  residuals  to  an  rms  of  ysec,  best  estimate 

of  the  frequency  drift  was  C.f-y  x  10’10  parts/day  which 
agrees  to  two  significant  places  with  the  irift  obtained  as 
a  by-product  of  the  Doppler  tracking  method. 
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Measures  of  the  accuracy  of  the  orbit  determination 
capability  of  the  passive  range  tracking  method  are  given  In 
Figs.  6.1  to  6.5.  The  H,  L,  C,  T  and  6  represent  the  correc¬ 
tions  as  a  function  of  time  that  must  be  added  to  the  ephem- 
erls  obtained  from  the  passive  range  data  to  obtain  the 
Doppler  determined  ephemerls.  The  K  refers  to  the  altitude 
difference,  the  L  to  the  along-track  difference,  the  C  to 
the  out-of-plane  difference  in  the  direction  of  the  orbital 
angular  momentum  and  the  T  represents  the  magnitude  of  the 
total  difference  which  is  equal  to  (H2  +  L2  +  C2)*'2.  The 
6  represents  the  angular  difference  In  the  two  ephemerides 
as  viewed  from  the  center  of  the  earth. 

The  total  difference,  T  ,  has  a  maximum  value  of 
59*0  nm  (109.2  km)  and  an  rms  of  30«9  nm  (57*3  km).  The 
angular  difference,  6  ,  has  a  maximum  of  2.7  milliradians 
and  an  rms  of  1.4  milliradians.  The  significant  difference 
in  the  two  ephemerides  Is  in  the  along-track  coordinate,  L. 
Since  this  difference  is  principally  a  linear  function  of 
time,  it  arises  from  a  difference  in  the  periods  (or  semi¬ 
major  axes)  in  the  two  ephemerides. 

Since  the  Doppler  determined  ephemeris  has  been 
Judged  in  [4]  to  be  accurate  to  better  than  100  nm  rms, the 
ephemeris  determined  from  the  passive  range  data  is  probably 
no  worse  than  [(3O.9)2  +  (100)2]V2  or  105  nm  rms.  This 
corresponds  to  an  angular  deviation  of  4,9  milliradians  rms. 
While  these  results  may  appear  to  be  rather  modest,  when 
compared  to  what  is  realizable  for  low  altitude  satellites, 
it  should  be  recalled  that  they  are  consistent  with  the 
original  objective  of  100  nm.  By  making  modest  improvements 
in  the  implementation  of  the  passive  range  tracking  system, 
the  technique  would  be  capable  of  providing  ephemerides  with 
a  substantial  improvement  in  accuracy.  Specifically,  it 
would  be  necessary  to  improve  the  synchronization  of  the 
clocks  at  each  of  the  ground  stations  and  to  utilize  either 
a  higher  frequency  or  a  dual  frequency  signal  to  minimize  the 
influence  of  the  ionosphere  on  the  time  of  propagation. 
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